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ABSTRACT 

Low-mass X-ray binaries (LMXBs) form efficiently in globular clusters (GCs). By combining Chandra 
X-ray Observatory and Hubble Space Telescope-Advanced Camera for Surveys observations of early-type 
galaxies, we probe the LMXB-GC connection using the most accurate identification of LMXBs and GCs to 
date. We explore the optical properties of 270 GCs with LMXBs and 6,488 GCs without detectable X-ray 
emission from a sample of eleven massive early-type galaxies in the Virgo cluster. Globular clusters that are 
more massive, are redder, and have smaller radii are more likely to contain LMXBs. Using known structural 
scaling relations for GCs, the latter implies that denser GCs are more likely to hold LMXBs. Unlike Galactic 
GCs, a large number of GCs with LMXBs have half-mass relaxation times > 2.5 Gyr; GCs do not need to 
survive for more than five relaxation timescales to produce LMXBs. By fitting the dependence of the expected 
number of LMXBs per GC, A;, on the GC mass M, color (g-z), and half-mass radius r/, coi, we find that 
A, cx M' '"**" "** 10''-'>^oa(s-z> r^2-2±0-2. This rules out that the number of LMXBs per GC is linearly proportional 
to GC mass (99.89% confidence limit). We derive an expression to estimate the number of multiple LMXB 
sources in GCs and predict that most GCs with high X-ray luminosities contain a single LMXB. The detailed 
dependence of A, on GC properties appears mainly due to a dependence on a combination of mass and radius, 
and a dependence on color, that are essentially equivalent to a dependence on the encounter rate F/, and the 



metallicity Z, A, cx F 



>0.82±0.05 70.39±0.{)7 



Our analysis provides strong evidence that dynamical formation and 



metallicity play the primary roles in determining the presence of an LMXB in extragalactic GCs. The shallower 
than linear dependence of GC sources requires an explanation by theories of dynamical binary formation; 
however, we note that our use of F/, as a proxy for the encounter rate, particularly if core-collapsed extragalactic 
GCs preferentially contain LMXBs, needs further testing in nearer galaxies. A metallicity-dependent variation 
in the number of neutron stars and black holes per unit GC mass, effects from irradiation induced winds, 
or suppression of magnetic braking in metal-poor stars may all be consistent with our derived abundance 
dependence; all three scenarios require further development. 

Subject headings: binaries: close — galaxies: elliptical and lenticular, cD — galaxies: star clusters — globular 
clusters: general — X-rays: binaries — X-rays: galaxies 



1. INTRODUCTION 

A low-mass X-ray binary (LMXB) consists of a compact 
stellar remnant, a neutron star (NS) or black hole (BH), that 
accretes the stellar envelope of a low-mass (M < 2Mq) stel- 
lar companion. Roche-lobe overflow is the major method of 
mass transfer in these systems, which in the Milk y Way have 
orbital periods of 0.19-398hr jWhite et al.lll"995D . When ac- 
tively accreting, LMXBs are bright X-ray sources 10^^- 
10^^ ergs"'; hereafter, the use of "LMXBs" refers only to 
active LMXBs). The Milky Way (L, - lO'^L©) contains 
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w 150 active LMXBs (iLiu et al.l2001h . Most of fliese LMXBs 
are in the field of the Milky Way and are likely to have 
been the result of the stellar ev olution of a primordial binary 
dVerbunt & van den H euvelT995). This primordial formation 
of an LMXB cannot be particularly efficient as the binary 
must survive the supernova that forms the NS/BH, yet remain 
in a tight enough orbit that Roche-lobe overflow eventually 
occurs. 

On the other hand, gravitational interactions in dense stel- 
lar systems, particularly globular clusters (GCs), could lead to 
more efficient formation of binaries in general, and LMXBs 
in particular. Over the history of X-ray astronom y, 13 ac- 
tive G alactic LM XBs hav e been detected in 12 GCs (iLiu et al.l 
'2001; White & Angelini 2001). In addition to active LMXBs, 
Galactic G Cs may contain seve n times as many LMXBs in 
quiescence (lHeinkeetal.ll2003h . Since GCs account for ~ 
0.1% of the Galactic fight, but 10% of the active Galactic 
LMXBs, formation of LMXBs in GCs must be hundreds of 
times more efficient than in the field (e.g., Clark' 1975; Kat^ 
Il975h . The increased efficiency may be due to a combination 
of two effects: (1) LMXB formation through tidal capture or 
exchange interactions after the supernova of the NS/BH; and 
(2) tighter binary orbits (binary hardening) in GCs due to ex- 
change interactions. The optical properties of GCs containing 
LMXBs provide a platform to test the dynamical formation of 
LMXBs in GCs. Furthermore, the metallicities of GCs can be 
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TABLE 1 
Galaxy Properties of Sample 



TABLE 2 

Properties of Chandra Observations 





Type 


''eff 


f SBF 


Mb 




Galaxy 


(RC3) 


(arcsec) 


(Mpc) 


(mag) 


Other Names 


(1) 


(2) 


(3) 


(4) 


(5) 




NGC 4365 


E3 


49.8 


23.2 


-21.46 


VCC731 


NGC 4374 


El 


51.0 


18.7 


-21.45 


M84, VCC763 


NGC 4382 


SO 


54.6 


17.9 


-21.46 


M85, VCC798 


NGC 4406 


E3 


104.0 


18.0 


-21.56 


M86, VCC881 


NGC 4472 


E2 


104.0 


16.8 


-21.98 


M49, VCC1226 


NGC 4486 


EO 


94.9 


15.7 


-21.70 


M87, VCC1316 


NGC 4526 


SO 


44.4 


16.1 


-20.60 


VCC1535 


NGC 4552 


EO 


29.3 


16.1 


-20.55 


M89, VCC1632 


NGC 4621 


E5 


40.5 


15.0 


-20.37 


M59, VCC1903 


NGC 4649 


E2 


68.7 


16.4 


-21.39 


M60, VCC1978 


NGC 4697 


E6 


72.0 


11.3 


-20.29 





readily estimated through their optical color, allowing one to 
probe for metallicity dependence in the evolution or forma- 
tion processes of LMXBs in GCs. Althoug h some patterns 
can be seen in the lirnited G alactic d ata (e .g. jBellazzini et alJ 
ll995l:lBregman et alj|2006l hereafter lBOel) . a larger sample is 
necessary to probe LMXB formation and evolution in GCs. 

With the launch of the Chandra X-ray Observatory iCXO), 
it is now possible to examine large samples of LMXBs. In 
particular. X-ray observations of nearby early-type galax- 
ies can identify tens to hundreds of bright X -ray sources, 
most of which are likely to be LMXBs (e.g., [Sarazin et al.' 
[2000 ). In these galaxies, - 20%-7 0% of the LMXBs appear 
to be associated wifli GCs (e.g.. iSarazin et al. 2000, 2001; 
lAngeUni etal]|2001l: iKundu et al.ll2002l: [Jordan et al. 2004b) ; 



early-type galaxies can provide the large number of GCs con- 
taining LMXBs necessary to probe LMXB formation and 
evolution in GCs. Early samples have shown that brighter 
and redder GCs are more likel y to c ontain LMX Bs (e.g . 
IKundu et allllool ISarazin et aLlBOOl iJordanet al.ll2004bh ; 
however, the use of resolution limited ground-based data or 
field-of-view (FOV) limited space-based data suppressed the 
numbers of GCs that could clearly be identified. 

By observing nearby early-type galaxies with the Wide 
Field Channel of the Hubble Space Telescope Advanced Cam- 
era for Surveys (HST-ACSfWFC), which has a 3f2 x 3.'2 FOV 
and can resolve GCs in nearby galaxies, one can construct 
accurate and comprehensive lists of GCs. In particular, the 
measurement of GC size allows better discrimination against 
unresolved background sources with GC colors, while pro- 
viding a GC property that is crucially important for probing 
the structure of GCs that contain LMXBs. The ACS Virgo 
Cluster Survey (ACS VCS) has ob served the centers of 100 
early-type galaxies (Cote et al. 2004) and provides excellent 
GC data. Some of its initial w ork explored the G C-LMXB 
connection in NGC 4486 (M87. lJordan et al.ll2004bl hereafter 
lACSVCSSh . In this work, we expand upon this analysis ex- 
ploring the GC-LMXB connection in 10 ACSVCS galaxies 
and a more nearby galaxy (NGC 4697) that may not belong 
to Virgo, but was observed following the ACSVCS observa- 
tional setup. We identify 270 GCs that contain LMXBs; this is 
one of the largest such samples to date, and the largest sample 
that includes information about GC size. 

2. SAMPLE 

2.L Galaxy Sample 

We selected our sample from early-type Virgo galaxies 
with detections of LMXBs using CXO and GCs using HST- 
ACS. Since the most massive galaxies are the most likely to 









'obs 


^exp 


Galaxy 


OBSID 


Detector" 


(ks) 


(ks) 


(1) 


(2) 


(3) 


(4) 


(5) 




2015 


API's S^/F 


40 4 


40 4 




0803 


rt-V^lkj O-J/ V 1 


28.5 


28 4 


NGC 4382 


2016 


ACTS-S3/F 


39.7 


39.7 




0318 




14 6 


11.4 




0963 




14 8 


12 


NGC 4472 


0321 


ACTS-S3WF 


39.6 


34.5 




0322 


AnS-TA/F 


10.4 


7 2 


NGC 4486 


0352 


ACIS-S3/G 


37.7 


33.6 




2707 


ACIS-S3/F 


98.7 


82.2 




3717 


ACIS-S3/F 


20.6 


9.2 


NGC 4526 


3925 


ACIS-S3A'F 


43.5 


41.5 


NGC 4552 


2072 


ACIS-S3A'F 


54.4 


54.4 


NGC 4621 


2068 


ACIS-S3/F 


24.8 


24.8 


NGC 4649 


0785 


ACIS-S3/VF 


36.9 


17.0 


NGC 4697 


0784 


ACIS-S3/F 


39.3 


37.2 




4727 


ACIS-S3/VF 


39.9 


39.9 




4727 


ACIS-S3/VF 


35.6 


35.6 




4729 


ACIS-S3AVF 


38.1 


32.0 




4730 


ACIS-S3A^F 


40.0 


40.0 



have large populations of both LMXBs and GCs, we con- 
centrated on galaxies with Mb < -20.25. Our sample in- 
cludes the ten brightest ACSVCS galaxies (NGC 4472, 4486, 
4406, 4365, 4382, 4374, 4649, 4526, 4552, and 4621) and 
NGC 4697, a similarly bright galaxy, which we observed 
with nearly the same setup as the ACSVCS. Table [T] sum- 
marizes some properties of the galaxies in our sample. The 
first three columns list the name, Hubbl e type, and effec- 
tive radius reft from ide Vaucouleurs et alJ (Il992h . The next 
columns list the (polynomial calibrated) surface brightness 
fluctuati on distan ce, Dsbf, determined from the HST obser- 
vations (iMei et al . 2006) and the absolute total B-magnitude, 
Mb (combining br from HyperLeda** and Z)sbf)- 

2.2. X-ray Analysis 

These galaxies all had archival or proprietary data from 
CXO; however, the observational setup varied widely between 
galaxies. Table |2] summarizes the properties of the CXO ob- 
servations, listing the galaxy name, CXO observation number 
(OBSID), detector with data mode (Faint, Graded, or Very- 
Faint), and the exposure times before (fobs) and after remov- 
ing periods of high backgrounds (fexp)- Since the ACIS-I chips 
are less sensitive than the ACIS-S3 chip, we only include data 
from ACIS-I when its fobs is > 25% the fobs of the ACIS- 
S3. Since the HST-ACS FOV is smaller than the CXO FOV, 
we hereafter only discuss X-ray sources within the CXO/HST 
overlapping region. 

All CXO observations were analyzed under CI AO 3.1'' with 
CALDB 2.28 and NASA's FTOOLS 5.3'". For Observa- 
tion 0784, the focal plane temperature was -110° C, while 
all other observations were taken at -120° C. When Very- 
Faint mode telemetry was available, the observations were 
cleaned using the extra data available in this mode to reduce 
the background level. All other observations were reduced in 
Faint mode, including OBSID 0352, which was telemetered in 
"Graded" mode. Known aspect offsets were applied for each 
observation. Our analysis includes only events with ASCA 
grades of 0, 2, 3, 4, and 6. Photon energies were determined 



»See 
9 See 



http : / / leda . univ-lyonl . f r /| 



http: //asc. harvard. edu/ciao/| 



' See |http : //heasarc . gsf c . nasa . gov/docs/ software/ lheasoft/| 
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using the gain files acisD2000-01-29gain_ctiN0001.fits, ex- 
cept for OBSIDs 0352 (acisD2000-07-06gainN0003.fits) and 
0784 (acisD1999-09-16gainN0005.fits). When appropriate, 
we corrected for time dependence of the gain and the charge- 
transfer inefficiency. All observations were corrected for 
quantum efficiency (QE) degradation and had exposure maps 
determined at 750 eV. We excluded bad pixels, bad columns, 
and columns adjacent to bad columns or chip node bound- 
aries. 

The use of local backgrounds in point source analysis mit- 
igates the effect of high background periods ("background 
flares"), which especially affect the backside-illuminated S 1 
and S3 chips". We avoided the periods with the most ex- 
treme flares by only including times where the background 
rate was below three times the expected rate from the blank- 
sky background fields in the CALDB. When available, we 
used the SI chip to measure the background rate, exclud- 
ing point sources, and compared to the blank-sky background 
count rates in CALDB in the 2.6-6. OkeV band. If the SI 
chip was not available, we used the S3 chip to determine 
count rates in the 2. 6-7. OkeV band. For Observation 0784, 
we compared count rates in the 0.3-10.0 keV band to Maxim 
Markevitch's aciss_B_7_bg_evt_271 103. fits blank-sky back- 
ground'^. When we binned the observations regularly in time 
to determine the count rate, small time intervals at the edges 
of the existing good-time intervals could not be included in a 
bin. This led to only a small amount of time that was lost in 
any observation. We list the exposure before and after flare 
removal in Table |2] 

To identify the discrete X-ray source population, we applied 
a wavelet detection algorithm (the CIAO wavdetect program) 
with scales increasing by a factor of \/2 from 1 to 32 pixels. 
We adopted a source detection threshold of 10"^ (< 1 false 
source per chip) for both the ACIS S3 and I chips. Source de- 
tection excluded regions with an exposure of less than 10% of 
that of the observation (e.g., regions at the edges of chips). To 
maximize the signal-to-noise (S/N), we analyzed the wavelet 
detection results from combined observations when available. 
We detected 708 X-ray sources that were also in the FOV 
of the HST-ACS observations. There are regions of complex 
gas emission at the center of NGC 4486. Although detected 
sources in these region may be LMXBs, they may also be 
compact regions of enhanced g as emission. To avoid mistaken 
identification of point sources, IACSVCS3I excluded aU X -ray 
sources in these regions. We follow this procedure, excluding 
33 X-ray sources. 

We used the coordinate list generated by WAVDETECT in 
ACIS Extract 3.34'^^ to refine the source positions and deter- 
mine source extraction regions. This was accomplished by 
determining the mean positions of events in source extrac- 
tion regions consistent with the X-ray point spread functions 
(PSFs) at the positions of the sources. For most sources the 
median photon energy was 0.6-2.6keV and we required the 
extraction region to encircle 90% of the X-ray PSF at 1 .5 ke V. 
We determined the PSF at either 0.3keV or 4.5 keV for the 
few sources whose median photon energy was either softer or 
harder. In a few sources where the 90% PSF extraction re- 
gions overlapped, we used a lower percentage of the PSF to 
define the extraction region. 

When there were multiple observations of a galaxy, we 

See http : //cxc . harvard . edu/ cont rib /maxim/ acisbg/| 
See footnote ml 

See |http : / /www, astro .psu ■ edu/xray/docs/TARA/ae_users_ 



matched the positions of detected sources in individual ob- 
servations to those present in the observation with the largest 
exposure. Absolute astrometric corrections were applied for 
all sources in all galaxies using a two-dimensional cross- 
corre lation techn ique to match sources from the Tycho-2 Cat- 
alog (lH0g et all 12000). 2MASS Point Source a nd Extended 
Sourc e Catalogs''*, and the USNO-B Catalog dMonet et alj 
l2003h . An astrometric correction (separately in RA and Dec.) 
was applied if the offset was statistically significant. 

We used a local background from an annular region whose 
area was approximately three times that of each source's ex- 
traction region; these local backgrounds include the diffuse 
emission from the host galaxy. To insure that the local back- 
ground regions did not contain any significant number of 
source counts, the inner radius of the background region was 
taken to be the radius encircling 97% of the PSF. In cases 
where background regions overlapped or fell along node/chip 
boundaries, we slightly altered these overlapping regions, pre- 
serving the ratio of source to background areas and ensur- 
ing that the source region and background region had similar 
mean exposures. For each of the sources, the observed net 
count rates, their errors, and S/N were calculated in the 0.3- 
6. OkeV band by stacking the observations for each galaxy, 
correcting for background photons, and dividing by the sum 
of the mean exposure over each source region. In NGC 4374 
and NGC 4486, a few of the sources detected with WAVDE- 
TECT had negative net count rates when determined in this 
way (two sources in NGC 4374, five sources in NGC 4486); 
we excluded these sources from further discussion. 

In order to convert count rates into energy fluxes, we fit 
a single emission model to the spectra of all of the LMXBs 
in all of the galaxies. For each observation of each galaxy, 
we extracted the cumulative spectra of essentially all of the 
LMXBs. We did not fit the spectra of LMXBs whose lu- 
minosities were determined at the <3a level or X-ray point 
sources that could be associated with both a GC and a non- 
GC optical source. To avoid contamination from background 
AGNs and foreground stars, we excluded all detected X- 
ray sources that were located within 1" of a non-GC optical 
source in the HST images. The source and background re- 
gions of candidate LMXBs were used to extract the spectra. 
We binned the spectra, requiring at least 25 total counts per 
bin, and only considered bins completely in the 0.5-10.0 keV 
band. We simultaneously fit the spectra of all of the sources in 
all of the observations of all of the galaxies to a single emis- 
sion model, which was taken to be either a power-law or ther- 
mal bremsstrahlung. However, we accounted for the differing 
Galactic absorbing columns (Nh) to the different galaxies, us- 
ing the Tuebingen-Bo ulder absorption (TB ABS) model assum- 
ing abundances from IWilms et al.l (|2000|) and photoelectric 
absorption cross-sections from lVerner et alJ (Il996). Note that 
the response files for each separate observation and source 
include the varying effects of absorption by the contaminant 
that produces the QE degradation in the ACIS detectors. The 
best-fit emission model for all of the sources was found to be 
a bremsstrahlung model with kT = 9.08;!;Q74keV. This model, 
combined with the individual values of Nh and distance for 
each galaxy and the individual response files for each observa- 
tion, were used to convert count rates into unabsorbed X-ray 
luminosities Lx in the 0.3-10. OkeV band. 

Due mainly to varying exposure times and distances, the 

'** When a source appeared in both 2MASS catalogs, the Point Source 
g^leg Bml'^ns were used. 
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TABLE 3 
GC-LMXB Matches by Galaxies 



Galaxy 


Nx 






fe-z)div 


^BlueGCs 


^RcdGCs 


^GCswLMXBs -'VbIucGCswLMXBs A'RedGCswLMXBs M'alseMatches 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) (9) (10) (11) 



Detected: All Sources With Positive Luminosity 



JNCjC 43 dD 


85 


0. 10 


yoo 


1.160 


418 [3] 


488 [7] 


43 (5) 


11 (0) 


29 (2) 


4.3d 


M/"' All A 

JNUL. 43 /4 


49 


0.05 


506 


1.200 


310 [0] 


196 [0] 


18(0) 


9 (0) 


9 (0) 




NUL. 43 o2 


37 


0.52 


505 


1.230 


320 [1] 


185 [1] 


11 (1) 


7 (0) 


3 (0) 


1 Ar\ 
i.4U 


JNLtL. 44Uo 


28 


0.09 


jo / 


1.220 


260 [0] 


107 [0] 


6(0) 


4 (0) 


2 (0) 


0.75 


JNLtL. 44 iZ 


86 


0.57 


764 


1.180 


306 [0] 


458 [4] 


40(2) 


4 (0) 


36 (2) 


2.o3 


Mr^r^ A A OA 
JNLtL. 44so 


84 


0. 14 


1639 


1.190 


690 [4] 


949 [10] 


49 (6) 


12 (0) 


34 (3) 


o.yD 


JNLiL. 4jzo 


34 


0.31 


244 


1.170 


121 [0] 


123 [0] 


7(0) 


3 (0) 


4 (0) 


0.98 


Mr^ /I ^ ^ T 
iNLiL. 4j JZ 


lo 




A^^ 


1.190 


217 [1] 


238 [1] 


31(1) 


o (U) 


22 (U) 


1 An 


Mnr^ /I AT 1 








1.110 


121 [0] 


185 [2] 


17(1) 


I (U) 




1 OA 




60 


65 


806 


1.200 


335 [0] 


471 [0] 


32(0) 


^ {(W 


zy 


1 






u.uo 




1.127 


127 [2] 


171 [2] 


34 (2) 




LI 


1 
Z. i J 










SNR: All Sources with 


Luminosities Determined at the 3 










29 


2 26 


906 


1.160 


418 [1] 


488 [5] 


17(3) 




1 s (1\ 


1 25 


NHP 4.'^74. 


16 


2.18 


506 


1.200 


310 [0] 


196 [0] 


5(0) 






0.62 


i>l V IV. JOZ 


21 


1 64 


SOS 


1.230 


320 [0] 


185 [0] 


5(0) 


d ((W 


1 

1 


u.oo 




19 


2.65 


367 


1.220 


260 [0] 


107 [0] 


6(0) 


d ((Yi 


z yyj) 


44 


NHP 4.479 


54 


1 30 


764 


1.180 


306 [0] 


458 [2] 


26(1) 


J [\J) 


zj (^i; 


1 72 




61 


1 30 




1.190 


690 [3] 


949 [5] 


40(3) 


Q ((W 


zy yv) 






12 


1 30 


244 


1.170 


121 [0] 


123 [0] 


3(0) 




z yyj) 


0.33 


NGC 4552 


49 


0.80 


456 


1.190 


217 [1] 


239 [1] 


25(1) 


6 (Q) 


18 fO^ 


1.53 


NGC4621 


16 


1.58 


306 


1.110 


121 [0] 


185 [0] 


4(0) 


0(0) 


4(0) 


0.55 


NGC 4649 


21 


2.85 


807 


1.200 


335 [0] 


472 [0] 


13(0) 


2(0) 


11 (0) 


0.43 


NGC 4697 


66 


0.15 


298 


1.127 


127 [2] 


171 [2] 


26 (2) 


7(1) 


19(1) 


1.75 










Complete: 


All Sources with Luminosities > 3.2 X lO'** 


ergs ' 






NGC 4365 


19 


3.20 


907 


1.160 


418 [1] 


489 [5] 


10 (3) 


1 (0) 


8(2) 


0.93 


NGC 4374 


12 


3.20 


506 


1.200 


310 [0] 


196 [0] 


4(0) 


2(0) 


2(0) 


0.45 


NGC 4382 


9 


3.20 


505 


1.230 


320 [0] 


185 [0] 


1(0) 


1(0) 


0(0) 


0.43 


NGC 4406 


16 


3.20 


367 


1.220 


260 [0] 


107 [0] 


5(0) 


3(0) 


2(0) 


0.38 


NGC 4472 


19 


3.20 


764 


1.180 


306 [0] 


458 [2] 


9(1) 


0(0) 


9(1) 


0.62 


NGC 4486 


27 


3.20 


1639 


1.190 


690 [3] 


949 [5] 


14(3) 


2(0) 


10(1) 


2.26 


NGC 4526 


4 


3.20 


244 


1.170 


121 [0] 


123 [0] 


1(0) 


1(0) 


0(0) 


0.11 


NGC 4552 


17 


3.20 


456 


1.190 


217 [1] 


239 [1] 


10(1) 


3(0) 


6(0) 


0.45 


NGC 4621 


7 


3.20 


308 


1.110 


121 [0] 


187 [0] 


2(0) 


0(0) 


2(0) 


0.23 


NGC 4649 


20 


3.20 


807 


1.200 


335 [0] 


472 [0] 


10(0) 


2(0) 


8(0) 


0.54 


NGC 4697 


3 


3.20 


298 


1.127 


127 [0] 


171 [0] 


3(0) 


0(0) 


3(0) 


0.00 



Units are 10"^^ ergs ^ in 0.3-10 keV band.!" If an X -ray source was within l"of GCs and non-GC sources, it is unclear if the X-ray source is 
matched to a GC; the possibly matching GCs have been removed from this sample. Numbers in brackets indicate the number of GCs where (only) 
multiple GCs are within l"of an X-ray source.*^ The number of RGC-LMXBs and BGC-LMXBs will not add up to the number of GC-LMXB s if 
an X-ray source is within 1 "of both an RGC and a BGC. Numbers in parenthesis indicate the number of X-ray sources where (only) multiple GCs 
are within l"of an X-ray source. 



observations of different galaxies have differing limiting sen- 
sitivities. For each galaxy, we summarize the number of X-ray 
detections and the minimum X-ray luminosity for three sam- 
ple definitions in columns 2 and 3 of Table [3] The three sam- 
ples are: all sources with a positive luminosity (Detected sam- 
ple); all sources with luminosities determined at the >3it level 
(Signal-to-Noise-Ratio [SNR] sample); and all sources with 
Lx > 3.2 X 10-"* ergs"' (Complete sample). All figures in this 
paper display the Detected sample. Note that the SNR sample 
does not include sources that were brighter than the reported 
minimum luminosity but not 3(t significant. Since a source 
with 20 net counts would be detected at high completeness for 
most of the galaxies, we defined the Complete sample using 
the luminosity that a source in NGC 4649 (the galaxy with the 
highest minimum luminosity) with 20 net counts would have, 
= 3.2 X 10-"^ ergs"'. We note that detections at the cen- 
ters of the X-ray brightest galaxies (e.g., NGC 4649) may still 
be incomplete due to the presence of bright gaseous emission; 
however, we have chosen this luminosity as a compromise be- 
tween completeness and a reasonable luminosity lower limit. 

2.3. Optical Analysis 



All of the galaxies in our sample, except NGC 4697, were 
observed as part of the ACSVCS, which acquired two 360s 
exposures in the F475W band (g-band), two 560 s expo- 
sures in the F850LP band, and one 90s F850LP exposure 
(z-band). D a ta redu ctions were carried out as described in 
[Jordan et al] (l2004al) . NGC 4697 was observed separately 
from the ACSVCS, but in a similar manner; its two F475W 
exposures were 15 s longer. We excluded optical sources de- 
tected in the central few arcseconds of som e galaxies with 
dusty cores (NGC 4374 and NGC 4526 in iFerrarese et"an 
l2006l and NGC 4697) and in the region of NGC 4649 that 
overlaps with its nearby spiral neighbor, NGC 4647. 

2.3.1. Observed GC Properties 

In our sample, over 10,000 optical sources were character- 
ized by their magnitudes, g and z, half-light radii, r/,, and 
po sitions as determin ed by KINGPHOT (see the Appendix 
in I Jordan et al.ll2005l where this code is described and sim- 
ulations are used to illustrate its performance). We note 
that ~ 80% of the GCs in our sample have r/, measured at 
SNR > 5. All magnitudes were converted to absolute mag- 
nitudes (Mg and Mr) using the surface brightness fluctuatio n 
distances in Table [T] As described in [Jordan et al.l (l2004al) . 
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Fig. 1. — GC magnitudes (M-) versus GC colors (g — z), with integrated liistograms of the properties above and to the right. GCs unmatched to LMXBs are 
indicated by small black dots and unfilled black histograms (scaled down by a factor of 10). Blue GCs with LMXBs are indicated by filled blue squares and 
histograms. Red GCs with LMXBs are indicated by filled red circles and histograms. The histograms of the GCs with LMXBs are stacked on each other. GCs 
that are redder and brighter are more likely to contain LMXBs. 



we first use magnitude, color and size criteria to select an ini- 
tial set of GC candidates. Then we use a statistical clustering 
method, described in detail in another paper in the ACS VCS 
series (Jordan et al. 2007, in preparation), which assigns to 
each source in the field of view of each galaxy a probabil- 
ity Pgc that the source is a GC. Our samples of GC candi- 
dates are then constructed by selecting all sources that have 
Pgc > 0.5. The res ults of our classification method are illus- 
trated in Figure 1 o fPeng et al.U200 6a). This selection returns 
7,084 likely GCs in our sample galaxies. We list the num- 
ber of GC candidates detected in each galaxy (excluding GCs 
whose matching X-ray source may also match to a non-GC 
optical source and all optical sources in the X-ray excluded 
regions of NGC 4486) in the fourth column of Table[3] 

Early-type galaxies often have a bimodal distribution of GC 
colors. For each galaxy, we used the division po int, ig-z)div, 
in the (g — z) color foUowing lPeng et aP (l2006al) to divide the 
GCs into blue-GCs and red-GCs; these division point ranged 
from (g-z)div = 1.11 to 1.23. We list the division points and 
the number of blue-GCs and red-GCs in columns five, six, and 
seven of Table [5] 

The positions of the GCs were used to determine the galac- 
tocentric distance, dec- These values were then scaled to the 
effective radii of the galaxies. We show scatter plots and in- 
tegrated histograms of the observed properties of the GCs in 
our sample in Figures[T]and|2l 



2.3.2. Derived GC Properties 

From the magnitude, colors, and radii, we can derive other 
parameters that might impact the formation and evolution of 
LMXBs in GCs. We calculate the GC mass (M) directly from 
the z-band magnitude, 

M=TAO-^^-^'-'^-'^-'^^Mq, (1) 

where is the z-band mass-to-light ratio and M-,q = 
4.512 is the absolute z-band magnitude of the Sun ob- 
tained from calcphot in the STSDAS Synphot IRAF pack- 
age. To estimate T,, we use vers ion 2.0 of the PEGASE code 
dFioc & Rocca- Volmerang3[T997ft. I n PEGASE, we used the 
stellar initial mass function of iKen nicutt (1983) to compute 
ig-z) colors and z-band luminosities for an assumed age 
r = 13 Gyr for several fixed values of [Fe/H]. Under these 
assumptions, is never more than 10% different from 1.5 
for [Fe/H] in the range -2 < [Fe/H] < 0, and thus z-band 
magnitudes are very good tracers of the mass. Comparably 
small ranges of T. result if younger GC ages are assumed. 
Given the average ig-z) of the GCs in any of our galaxies, 
we interpolate on these PEGASE model curves to estimate 
average z-band mass-to-light ratios. The values of T, depend 
slightly on the mean color of the GC systems, but the variation 
was less than 1 % across our sample galaxies when assuming 
that GCs are uniformly old. Thus, we chose to use a constant 
T, = I.45M0/LQ for all galaxies. 
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Fig. 2. — GC galactocentric distances (doc) versus observed GC half-light radii (r/,), with integrated histograms of the properties above and to the right. The 
symbols follow those of Figure[T] GCs that are smaller in extent are more likely to contain LMXBs. The galactocentric distance does not appear to affect whether 
a GC contains an LMXB. 



There is an observed correlation between r/, and (g-z), 
in the sense that metal-rich (red) GCs are found to be w 
17% smaller than their metal-poor (blue) counterparts (e.g., 
[Jordan etalJ[2()05h . This can be a consequence of mass- 
segregation combined with the metallicity dependence of stel- 
lar lifetimes under the assumption that the avera ge half-ma5 5 
radius of GCs does not depend on metallicity (lJordanll2004l) . 
The observations of the size difference between metal-rich 
and metal-poor GCs in V irgo are consistent with this explana- 
tion (iJordan et al.ll2005 l). Thus, it is possible that the smaller 
half-light radius for metal-rich GCs do not reflect a decrease 
in the corresponding half-mass radius r/, m- 

Since ri,_M is more physically relevant in discussion of 
Galactic dynamics than r/,, we calculated a "corrected" half- 
light radius as 

rH.cor = n , 10" '^[<«-^)-'-21 = r„,M , (2) 

following [Jordan et al.l (l2005l) . and assume in what follows 
that this is the GC half-mass radius. Note that while this 
definition takes care of a color dependence in the half-light 
radii that might not reflect a corresponding dependence in the 
half-mass radii, r/, cor can still differ from the actual half-mass 
radius by a constant factor (Figure 2 in Jordan 2004). We 
display flie r^ cor of the GCs in our sample in Figure |3] 

Since lB06l suggested that LMXBs occurred in GCs with 
ages greater than five times the relaxation time at the half- 
m ass radius, f/; . reiax, w e adopted the Harris, (1996) corrections 
to iDiorgovskil (1 19931) equation (11) and used our corrected 



half-mass radius estimate for calculating f/7.reiax.c 



^/j, relax, cor ~ 



2.055 X lO*- 



ln(0.4A^*) V ^< 



M 



1/2 



Ipc 



3/2 



yr, 

(3) 

where the average stellar mass, (m*), is taken to be \Mq and 
N^=M/{m^). B06 suggested that the relaxation time be com- 
bined with the number of stars in the globular cluster to give 
a "stellar interaction rate" of 



^/7, relax, C' 



■yr (xM ' r 



1/2,-3/2 



h.M 



(4) 



For a system in virial equilibrium, the relaxation time is 
approximately (0.1A^*/lnA^*)fcros.s, where fcross is flie stellar 
crossing time (roughly the orbital period of a star within the 
GC). Thus, S w 101n(A^*)/f(;ro,ss; except for a nearly constant 
factor, it is just the inverse of the orbital period. We will refer 
to S as the "stellar crossing rate." In Figure|4] we plot f/,.reiax.cor 
versus M- for the GCs in our sample, overlaying lines of con- 
stant S for comparison with B06, Figure 5. 

The dynamical formation model of LMXBs in GCs is the 
leading explanation for the larger efficiency of LMXB pro- 
duction in GCs compared to the fields of galaxies. In this 
model, the binary is either formed by tidal capture or through 
exchange interactions between a NS/BH and an existing bi- 
nary. In both cases, the encounter rate (F) is thought to depend 

on the properties of GCs at their cores; F cx f^J^r^, where po is 
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Fig. 3. — Estimated GC half-light radii (r;, cor) versus GC colors (g—z), with integrated histograms of the properties above and to the right. The symbols follow 
those of Figure[T] There is no color dependence of r/, cor. indicating we have successfully removed the color dependence of half-light radius r;,. GCs that are 
redder and smaller in extent are more likely to contain LMXBs. 

the central density, is the core radius, and the virial theorem 
has been assumed to c onnect core velocities to densities and 
radii (see, e.g.. I Verbuntii2003() . Although these parameters are 
not directly measurable at the distance of the Virgo cluster, 
we can create a proxy. 



r/,= 



/ M 



3/2 



Ipc 



-5/2 



(5) 



where we have assumed that the relation between core struc- 
tural parameters and the parameters at a half-mass radius (dic- 
tated by the concentration, c) is the same for all GCs. The 
latter assumption is not necessarily true, e.g. there is an ob- 
served trend in Galactic GCs for the concentration to increase 
with GC mass (McLaughlin 2000). In any case, we define 
F/, to be our tracer for encounter rates under the assumptions 
stated; our results on the dependence of the expected number 
of GC LMXBs on F/, can be easily compared with theoreti- 
cal predictions based on any alternate assumptions. We note 
that if GCs with higher concentrations are more likely to con- 
tain LMXBs, then the F/, we calculate will underpredict the 
encounter rate for GCs with LMXBs. We plot r/, cor versus 
M, for GCs in our sample, ov erlay ing lines of constant F/„ in 
F igure [5]for c omparison with lBOa, Figure 3. 

lACSVCSSl used measured King models concentrations c to 
estimate F directly. As the measured c are rather uncertain 
for most sources, we prefer to use F/, in this work. We note 



though that the results of IACSVCS3I are robust with respect 
to uncertainties in the measured c; indeed all conclusions re- 
main unchanged if a single c = 1.5 is assumed for all GCs. 
In other words, their conclusions remain the same when using 
F/, rather than F. 

2.4. Matching LMXBs And GCs 

To determine the relative astrometry between CXO and HST 
observations for each galaxy, we convolved the offset posi- 
tions in RA and Dec between all X-ray and optical positions 
with a a = 0."5 2-dimensional Gaussian, which approximates 
the PSF of CXO, to create a cross-correlated image. The max- 
imum cross-correlation was used to determine the astrometric 
offset. After correcting for this offset. X-ray sources within 
1" of optical sources were considered to be matched. We 
chose 1" as a balance between accurate identification of the 
X-ray source astrometry (the astrometric accuracy of CXO in 
the small count regime is ^ 0.3-0.5") and the increased num- 
ber of falsely identified sources with larger match radii. 

If an X-ray source was matched to both non-GC optical 
sources (e.g., background galaxies) and GCs, we excluded 
the X-ray source and GCs from our analysis. All remain- 
ing 288 X-ray detections in GCs were considered to be GC- 
LMXBs. For the 18 GC-LMXBs that were within 1" of 
multiple GCs, we could not determine which GC contained 
the LMXB; however, if the multiple GCs all belonged to the 



8 



SIVAKOFF FT AL. 



-— I ■ 1 ■ 1 1 1 ■ f 




-12 -10 -8 -6 
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displa yed in iBO& Figure 5. GCs with large relaxation timescales can host LMXBs, in contrast to the result found in the much smaller sample of Galactic GCs 
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Fig. 5. — Estimated GC size ifh,coT) versus GC magnitude (M-), where the symbols are the same as in Figure \T\ O verlaid are lines of constant encounter rate 
that decrease by 0.5 dex starting at logF;, = 9.0 (lower left). The Galactic version of this plot for V-band is displayed in lBOd Figure 3. GCs with larger encounter 
rates are more likely to contain LMXBs. 
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same sub-population (blue-GC vs. red-GC), we could count 
the GC-LMXB as belonging to that sub-population. We indi- 
cate the number of GCs with LMXBs, and separately list the 
numbers of LMXBs in blue-GCs and red-GCs, in the eighth, 
ninth, and tenth columns of Table [3] In parenthesis, we indi- 
cate the number of GC-LMXBs matched to multiple GCs. We 
also indicate the corresponding number of GCs in the brackets 
of columns four, six, and seven. 

We separated the GCs into two sub-populations, those that 
clearly contained an LMXB in our Detected sample (indicated 
by a subscripted X), and those that clearly did not (indicated 
by a subscripted nX). We found 270 GCs that contained an 
LMXB at our detection levels and 6,488 GCs that clearly did 
not. We display the properties of these GCs in Figures 
indicating blue-GCs with LMXBs using filled blue squares 
and red-GCs with LMXBs using filled red circles. Of the 270 
GCs clearly containing an LMXB, 160 are in the SNR sample, 
and 61 are in the Complete sample. 

We randomized the position angles of the GCs around 
the centers of their host galaxies to estimate the number of 
LMXBs falsely associated with GCs. We adopted two differ- 
ent methods: either we kept the galactocentric distance of the 
GC fixed, or we assumed that the GC lay on the same elliptical 
isophote on which it was originally detected. We reassessed 
the matches with the randomized GC catalogs, and used the 
number of resulting matches to predict the percentage of non- 
GC X-ray sources falsely matched to a GC, /?faise- Since there 
was typically little difference between the two methods, we 
averaged the two. Given the typical values of these num- 
bers (3.3-14.8%), it is unlikely that any of the GC-LMXBs 
matched to multiple GCs are not actually associated with a 
GC; however, some fraction of the GC-LMXBs matched to a 
single GC will be false matches. The predicted number of 
falsely matched sources is /7faise/(l "/'false) times the num- 
ber of observed non-GC X-ray sources. The term in the de- 
nominator accounts for the falsely matched sources removed 
from the non-GC X-ray sample. The expected number of 
false matches is displayed in the eleventh column of Table [3] 
The large GC density in NGC 4486 leads to a correspond- 
ingly larger number of false matches; however, the number 
of false matches still make up less than 15% of the identified 
GC-LMXBs. With the number of false GC-LMXBs and to- 
tal number of GCs, we estimate the expected number of false 
matches per GC, A /, for each galaxy, which range from about 
2 X 10"^^ to 7 X lQ-\ 

3. PROPERTIES OF GCS WITH AND WITHOUT LMXBS 

In what follows we compare the properties of GCs that har- 
bor an LMXB with those of GCs that do not harbor one. We 
probe in turn the distributions of luminosity (mass), (g - z) 
color (metallicity), size (r/,), galactocentric distance, half- 
mass relaxation time f/,.i-ei and dynamical rates S. By revealing 
which properties drive the presence of an LMXB in GCs - or 
which ones are irrelevant - this exercise can shed light on the 
formation process of LMXBs in GCs. 

We employ three methods to compare the properties of 
GCs with and without LMXBs (When comparing the proper- 
ties of the SNR sample and Complete sample to GCs with- 
out LMXBs, we did not add the GCs which had LMXBs 
that did not fall in the sample into the sample of GCs with- 
out LMXBs.) First, we use binned histograms of various 
properties to display the qualitative differences. These his- 
tograms are displayed in Figures [T]-0] Second, we calcu- 
late the median values of the GC properties and use the 



non-parametric Wilcoxon rank-su m test (fJwRS, equivalen t to 
the Mann- Whitney rank-sum test: lMann & Whitneylll947h to 
quantify the differences. We additionally use the Wilcoxon 
test to compare the properties of GCs with LMXBs in the 
Complete sample and GCs with LMXBs fainter than the 
Complete sample limit. Finally, we bin the GCs by the dif- 
ferent parameters available and look for non-uniform proba- 
bilities of GCs containing an LMXB. The results of the latter 
exercise are shown in Figure |6] The errors on the fractions in 
Figure|6]are calculated assuming Poisson statistics (i.e., strict 
la confidence intervals were calculated as opposed to using 
the a/^V approximation). The fractions in Figure |6] were cor- 
rected for the rate of false matches. A/. 

3.1. Luminosity and Mass 

Prior observations have revealed that more luminous GCs 
appear to preferen t ially h ost LMXBs (e.g., A ngelini et all 
2001; Kundu et all 120021: ISarazin et all 120031; Ijordinetail 
i2004b(). Figure [T] confirms that LMXBs are found more of- 
ten in brighter GCs. The median is -8.5 for GCs with- 
out LMXBs, and -9.9 for GCs with LMXBs (all samples), 
which correspond to Wilcoxon rank-sum differences of 15.6ct 
(Detected sample), 1 1.8cr (SNR sample), and 6.6(j (Complete 
sample). We find no significant difference (ctwrs = 0.9) in the 
optical luminosities of GCs in the Complete sample and of 
GCs with fainter LMXBs. 

The median mass of a GC with an LMXB is 3.6 times that 
of a GC without an LMXB. The upper left panel of Figure |6] 
suggests a power-law dependence of the probability of a GC 
containing an LMXB on the mass of the GC. 

3.2. Color 

Prior observations also have revealed t hat redder GCs ap- 
pear to preferentially host LMXBs (e.g.. Ku ndu et al.l 120031 ; 
Sarazin et al. 2003; Jordan et al. 2004b). Our larger sample 
clearly confirms this (see Figure[T]). The median (g-z) color 
is 1.20 for GCs without LMXBs, and 1.34 (Detected sample), 
1 .35 (SNR sample), and 1.38 (Complete sample) for GCs with 
LMXBs, which correspond to Wilcoxon rank-sum differences 
of 6.6(7, 5.8(7, and 3.7(7. There were no significant differences 
(cwRS = 1 -2) in the colors of GCs in the Complete sample and 
GCs with fainter LMXBs. 

Previous works have often compared the fraction of blue- 
GCs with LMXBs to the fraction of red-GCs with LMXBs, 
and found that red-GCs are ^ 3 times more likely to contain 
LMXBs that their blue counterparts. In Figure|7l we compare 
these fractions computed using the Detected sample. The cu- 
mulative sample suggests that red-GCs are 3.15 ± 0.54 times 
more likely to have LMXBs than blue-GCs; however, there is 
considerable scatter in this fraction between galaxies. 

The upper right panel of Figure |6] suggests that the proba- 
bility that a GC contains an LMXB depends exponentially on 
the color. Based on the observed color distribution bimodality, 
some GC formation scenarios posit that the red and blue GC 
subpopulations h ave different formation histories (see, e.g., 
IWest et alj |2004). If the presence of an LMXB in a GC was 
somehow mainly determined by this history, it is conceivable 
that the probability of a GC holding an LMXB was differ- 
ent for each subpopulation but independent of color within 
each subpopulation. If this dependence were due only to 
a difference between the red and blue GC populations, but 
were independent of the color of a GC within the two pop- 
ulations, this would appear as a step-function relation in the 
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Fig. 6. — The percentage of GCs containing an LMXB as afunction of GC magnitude (Mj), color (g-z), half-mass radius {ri,,coi), and galactocentric distance 
(doc)- GCs that are brighter, are redder, and have smaller half-mass radii are more likely to contain LMXBs. The uncertainties are at the la confidence level. 
The bins were set by requiring 27 GCs with LMXBs per bin. 




2 4 6 8 

Blue-GCs With on LMXB [%] 

Fig. 7 . — The percentage of red-GCs containing an LMXB versus the percentage of blue-GCs containing an LMXB for the Detected sample. The data points 
with la confidence limits are the different galaxies in our sample. The dotted Une indicates the cumulative relation that red-GCs contain as many LMXBs as 
blue-GCs; the data clearly do not support this relation. The soUd Une indicates the cumulative relation that red-GCs contain 3. 15 ± 0.54 times as many LMXBs 
as blue-GCs. The greyscale indicates the area covered by 1, 2, and 3a off of the best fit. 
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upper right panel of Figure |6] except for the third bluest color 
bin. This bin contains all but two of the different galactic 
division points between red-GCs and blue-GCs. As such, 
its value would be intermediate between the two efficiencies. 
Since both pictures, i.e., a continuous exponential dependence 
and a step-function, have two constraints, we can compare 
their ^ fits to test which interpretation is more consistent 
with the data. For these fits, we always excluded the third 
bluest color bin and tested the effect of excluding zero, one, or 
two of the reddest color bins. Since the step-function model 
always had a higher compared to the exponential func- 
tion model (x^(j.p/dof = 12.8/6, 12.8/5, and6. 2/4 compared 

to Xexp/dof =9.7/6,6.0/5, and 3. 7/4, respectively), we do not 
believe that the probability that a GC contains an LMXB de- 
pends on a difference between the red and blue GC popula- 
tions that is independent of the color of a GC within the two 
populations. Since color is roughly a logarithm of metallicity, 
we prefer interpreting the upper right panel of Figure |6] as an 
indication that the probability a GC contains an LMXB has a 
continuous power-law dependence on GC metallicity. 

3.3. Size 

Although there were indirect indications that the size 
of GCs af fect the formation or evolution of LMXBs in 
|ACSVCS3[ we present here the first explicit evidence of this. 
We find that LMXBs are found more often in GCs that have 
smaller half-light radii (Figure |2]i. Becaus e the hal f-light radii 
of GCs are uncorrelated with their mass dMcLaughhn 2000; 
[Jordan et al.ll200 5). this implies that LMXBs are found more 
often in GCs that are denser. The median r/, values are 2.6 pc 
for GCs without LMXBs, and2.3pc (Detected sample), 2.3pc 
(SNR sample), and 2.2 pc (Complete sample) for GCs with 
LMXBs, which correspond to Wilcoxon rank-sum differences 
of 8.1(7, 7.4(7, and 4.9(7. Although the differences in the me- 
dians are < 15%, the distributions are found to be signifi- 
cantly different for two reasons. First, with the number of 
GCs in our samples, the Wilcoxon test is sensitive to 10%, 
14%, and 22% differences in the medians, respectively, at the 
3(7 level. Second, the Wilcoxon test samples the entire dis- 
tribution, not just the medians. For example, the values of 
the half-light radii that contain 90% of the GCs are 4.3 pc for 
GCs without LMXBs, and 3.3 pc (Detected sample), 3.0pc 
(SNR sample), and 3.0pc (Complete sample) for GCs with 
LMXBs. Kolmogorov-Smirnov tests indicate that the proba- 
bility the distributions of half-light radii are drawn from the 
same sample are 5.7 x lO"'-' (Detected sample), 8.1 x 10~" 
(SNR sample), and 9.6 x 10"'' (Complete sample), consistent 
with the Wilcoxon tests. Once again, we found no significant 
differences (ctwrs = 1 -2) in the half-light radii of GCs in the 
"Complete sample" and GCs with fainter LMXBs. 

Given that redder GCs are more likely to contain LMXBs 
and have smaller half-light radii, we also compared the half- 
mass radii (see eq. ||2l) of GCs with and without LMXBs 
(Figure |3]l. Our adopted half-mass radii should be, by defini- 
tion, nearly insensitive to color variations. The median r/, coi 
of GCs with LMXBs were only 0.1 pc more than their half- 
light radii. The distributions of half-mass radii of GCs with 
and without LMXBs are different at the 5.9(7 (Detected sam- 
ple), 5.5(7 (SNR sample), and 3.7ct (Complete sample) levels 
according to the Wilcoxon rank-sum test. Thus, GCs that are 
smaller in extent remain more likely to contain LMXBs, even 
if the half-mass radii are used to measure the size. 

The lower left panel of Figure |6] suggests that the probabil- 



ity that a GC has an LMXB may decrease roughly as a power- 
law of the half-light radius. There is a similar, but slightly 
flatter dependence on half-mass radius. 

3.4. Galactocentric Distance 

iKim et akl (|2006 |) present a study of LMXBs in six early- 
type galaxies (four overlap with our study). Combining HST- 
WFPC2 and ground-based observations they identify 285 
LMXBs with GCs. They find that LMXBs are more likely 
(at the 6.6(7 significance level) to be found in GCs in the 
inner 0.5iieff compared to an annulus between one-half and 
one times the isophotal B = 25 mag arcsec"^ radius; 44 GC- 
LMXBs out of 1004 GCs are found in the inner region versus 
69 GC-LMXBs out of 2908 GCs in the outer). The suggest 
their results indicate that galactocentric distance plays a criti- 
cal role and that GCs may have more compact cores near the 
galactic center, and thus higher encounter rates. 

The effect of galactocentric distance on the probability that 
a GC has a LMXB is less clear in our data. For compari- 
son with their result, we examine the same inner annulus and 
1 .5(ieff <d < 3.Qdss, which corresponds to approximately the 
same outer annulus. We find 66 GC-LMXBs out of 121 1 GCs 
versus 44 GC-LMXBs out of 1203 GCs. Although the in- 
ner region is numerically more likely to have GCs containing 
LMXBs in our sample, this result is not statistically signifi- 
cant (1.7(7). Comparing the inner regions of the two samples, 
where their study relies heavily on HST identifications, are 
clearly consistent. There is a large discrepancy in the outer 
region. With ground-based (and to a much lesser extent HST- 
WFPC2 observations), it is more likely that contamination 
of GCs with unrelated objects occurs w hen com p ared w ith 
our HST-ACS observations. Although Kim et alJ (l2006h at- 
tempt to account for this contamination, underestimation of 
the contamination level would result in a larger falsely iden- 
tified GC population that would deflate their measured likeli- 
hood of finding a GC with an LMXB in the outer field. 

In our data, the values of dec are slightly smaller for GCs 
with LMXBs in the Complete sample as compared to GCs 
with LMXBs not in the Complete sample, but the difference 
is not very significant (ctwrs = 2.2). This could indicate that 
the values of dec are affected by incompleteness. Since the 
incompleteness is more important in the X-ray brighter, cen- 
tral regions of galaxies, and incompleteness affects the fainter 
LMXBs more, this effect is not unexpected. There is no sig- 
nificant difference ((7wrs = 1 -6) in the radial distributions be- 
tween the Detected sample of GCs with LMXBs and GCs 
without LMXBs; however, there are slightly significant dif- 
ferences when using the SNR sample ((7wrs = 3.1) and Com- 
plete sample ((7wrs = 2.8). The lower right panel of Figure|6] 
suggest that there may be a slight decrease with dec in the 
probability that a GC contains an LMXB. However, since red- 
GCs tend to be closer to the galaxy center and are more likely 
to contain LMXBs, it is possible that the lower galactocen- 
tric distances of GCs with LMXBs in the SNR and Complete 
sample can be completely accounted for without invoking an 
intrinsic difference in the probability a GC contains an LMXB 
on dec- Given the marginal, if any, strength of the dependence 
on lice, we do not consider it further in this paper. We have 
however checked that including a galactocentric distance ef- 
fect does not alter any of our scientific conclusions. 

3.5. Relaxation Time 

No Galactic GC with a rela xation timescale > 2.5 Gyr con- 
tains an active LMXB ( lB06h . Their Galactic sample is lim- 
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ited to 141 GCs, 12 of which contain LMXBs with Lx > 
10^^ ergs"'. Our sample, which does not probe as deeply 
in X-ray luminosity, includes 6,758 GCs, 270 of which con- 
tain LMXBs. In our much larger sample (Figure |4|i, about 
18% of the GCs with LMXBs have relaxation timescales 
larger than this. We used Monte Carlo simulations to esti- 
mate that measurement errors only bias this result by about 
3%. In other words, the relaxation timescale of about 15% 
of the GCs are intrinsically larger than 2.5 Gyr The relax- 
ation timescale of the GCs with LMXBs in our sample are 
larger than those GCs without LMXBs; the median f/i.reiax.cor is 
1 .OGyr for GCs without LMXBs, and 1 .5 Gyr (Detected sam- 
ple), 1.4Gyr (SNR sample), and 1.3 Gyr (Complete sample) 
for GCs with LMXBs. This difference is opposite to that seen 
in the Galaxy, where the median f^.ieiax are 0.7 Gyr for GCs 
with LMXBs and 1.3 Gyr for GCs without LMXBs (Harrii 
119961) . In our sample, the relaxation timescale distributions 
of GCs with LMXBs are different from GCs without LMXBs 
at the 6.8(7, 4.6cr, and 2.4ct level according to the Wilcoxon 
rank-sum test; the equivalent test on Galactic GCs indicates 
a difference of 2.0(7. We found no significant differences 
(fwRS = 1-3) between the half-mass relaxation timescales of 
GCs in the "Complete sample" and GCs with fainter LMXBs. 

3.6. Dynamical Rates 

There are two dynamical rates that could affect formation or 
evolution of LMXBs: the stellar crossing rate {S oc M'/^ r"^/^) 
and the encounter rate (F/, oc M^/^ r~^l^). Since GCs that 
are more massive and smaller in extent are more likely to 
contain LMXBs, it is clear that both of these rates would 
be higher in GCs with LMXBs compared to GCs without 
LMXBs (Fig. |5] and nil. The median values of 5 are 0.74kyr"' 
for GCs without LMXBs and 1.8kyr"' (Detected sample), 
2.0kyr~' (SNR sample), and 1.8kyr~' (Complete sample) for 
GCs with LMXBs, corresponding to \6.9a, 13.3(7, and 7.8ct 
Wilcoxon rank-sum differences. For F/,, the median values 
are 6.9 x 10^ for GCs without LMXBs and 7.3 x 10'' (De- 
tected sample), 7.9 x 10*" (SNR sample), and 6.4 x 10'' (Com- 
plete sample) for GCs with LMXBs, corresponding to 17.1(7, 
13.2ct, and 7.7(7 Wilcoxon rank-sum differences. We found 
no significant differences ((7wrs = 0.1) between either of the 
dynamical rates of GCs in the Complete sample and GCs with 
fainter LMXBs. 

4. MULTI- VARIABLE RELATION BETWEEN LMXBS AND GCS 
4.1. Technique 

Since several properties appear to affect the probability that 
a GC contains an LMXB, it is important to simultaneously 
account for these properties. Based on Figure |6] and the rel- 
atively small percentage of GCs falsely matched to LMXBs 
A/, we have assumed the expected (true) number of LMXBs 
in a GC (A,) has the following dependence on GC properties: 
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where we fit for the normalization (A) and indices ((3, 5, and 
e). Here, r can either be the measured half-light radius r;, 
or the corrected half-mass ra dius r/, cor- O ur approach differs 
from the fitting performed in lACSVCS3 l; we choose to fit the 
expected numbers of LMXBs per GC, as opposed to the prob- 
ability a GC contains an LMXB. The probability must saturate 
at unity, whereas the expected number can be unUmited; this 
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makes it easier to fit simple functions (e.g., power-laws) to the 
expected number and compare with Galactic results. The two 
approaches give the same results when A? is small. By includ- 
ing the normalization A we get a direct estimate for the num- 
ber of LMXBs in GCs and we are able to probe quantitatively 
the fraction of GC X-ray point sources that are expected to 
result from the integrated X-ray flux of multiple sources. An- 
other difference is that we include a term to account for GCs 
falsely matched to LMXBs, A = A/ + A,. 

The expected number of LMXBs in a GC can be converted 
to a probability that there are no LMXBs assuming a Poisson 
distribution, 

P„x = exp(-A), (7) 

and the probability that there is at least one LMXB, 

Px = l-exp(-A). (8) 

One can then vary the parameters to maximize the log likeli- 
hood. 
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where the products are taken over the lists of GCs with no 
LMXBs and GCs with LMXBs. Given that we did not find a 
difference in the masses, colors, or sizes of GCs with LMXBs 
in the Complete sample as compared to those GCs with fainter 
LMXBs, we first choose to use the better statistics provided 
by the Detected sample to fit the indices. Given the widely 
different luminosity limits in the Detected sample, we do not 
believe the normalization we fit is physically meaningful. To 
derive a physically meaningful normalization, we apply our 
derived exponents to the Complete sample. Since the log 
likelihood can be related to Ax , (-2Ai/) = Ax^), we use the 
change in log likelihood for one degree of freedom (dof) to 
determine one-dimensional fitting errors (1(7) on each vary- 
ing parameter. 

We display the results of our fits to the Detected sample in 
Table m where bracketed numbers indicate that an index was 
frozen to that number. Once a form for A, as a function of 
mass, color, and radius has been determined, we can use this 
expectation to predict the fraction of GCs having an LMXB 
as a function of each of the GC properties. These predictions 
can be compared to the plots of the variation of this proba- 
bility with individual properties shown in Figure |6l For each 
of the binned data points in each of the panels of this Figure, 
we apply the best-fit expectation calculated for each of the 
other properties of the GCs in the binned data point, and then 
compute the expected fraction integrated over that bin. This 
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is compared to the actual fraction for the bin, and the con- 
tribution to from that bin is determined. Since we do not 
consider the dependence of A, on galactocentric distance, we 
did not include that panel's contribution to x^. Since the x^ 
of that panel ranged from 4.2-1 1 .4 for 10 dof for all fits in Ta- 
bled this is justified. In this way, we can determine whether 
the best-fit expectation function fits the binned data in Fig- 
ure |6] adequately. The results, which we describe in detail in 
what follows, are shown in Figures[8]and|9] 

4.2. Effect of GC Mass and Color 

We first established a baseline assuming that the expected 
number of LMXBs per GC has no dependence on GC prop- 
erties (Tabled row 1); our binned x^ was 402.4 for 29 dof 
(3 X 10 data points minus one constraint for the normaliza- 
tion). We next individually varied the mass and color indices 
(Table |4] rows 2 and 3); mass clearly had the largest effect 
(-2A-0 = 290.2); however, its binned x^ was still as large as 
158.6 for 28 dof. If we included the effects of mass and color 
simultaneously (Table |4] row 4), we had a much improved fit 
from the mass-only fit (-2Aip = 37.6). We display this fit in 
Figure[8] which has a binned x^ of 106.4 for 27 dof. Although 
the effect of color on X, was included in this fit, the resulting 
expression does not provide a very good fit to the dependence 
of the fraction on color alone (Figure [8] upper right panel). 
The x^ is 26.5 for 10 bins in this figure. In general, the depen- 
dence of A, on mass and color does not adequately describe 
the observed number of LMXBs per GC as a function of the 
observed GC properties. Some dependence on another GC 
property is required. 

We note that based on known correlations of GC properties 
we do not expect a relation based on GC mass and color alone 
to be able to reproduce any dependence on GC size. Namely, 
the mean of GCs is observed to be independent of mass 
( lMcLaughlinll2000l: iJordan etal]|2005h and our r^ cor has no 
dependence on color. Thus, any dependence of A, on r/, ^of, 
as observed in the lower left panel of Figure |6] cannot be re- 
produced by a function which depends only on GC mass and 
color 

In the prior largest sample of GCs wi th LMXBs t o date 
(6 galaxies, 98 LMXBs in 2,276 GCs), ISmits et al.1 (|2006|) 
found P(LMXB) oc Mi'^^O'^Z^^^ioo^ since they used 
multiple color indices, they adopted linear color-metallicity 
relations based on Milky Way GCs, ignoring scatter in 
the relation. For comparison with their data, we adopt 
the linear color-metallicity relation presented in IACSVCS3I 
(g-z) ^ 0.38[FE/H]+ 1.62. In that case, we derive A oc 
^i.oo9_„(,5,20-29±o.io^ which is consistent with their results. 
However, as we point out in the next section, the effect of 
GC sizes must be taken into account. This effect alters the 
measured indices for mass and color, which will in turn alter 
interpretations based on those indices. 

4.3. Effect of GC Size 

There are two GC sizes we considered in our fits: the half- 
light radius (r = r/,), and the half-mass radius (r = r/, coi-)- Fit- 
ting half-light radius alone (Table |4] row 5) leads to a better 
fit than fitting color alone; however, given the correlation be- 
tween color and half-light radius, it is not clear if this is due 
in part to the dependence on color or a real dependence on 
the cluster size. On the other hand, only fitting the half-mass 
radius (Table ID row 7) is nearly equivalent to only fitting the 
color (row 3). This is consistent with the picture developed 
from the Wilcoxon tests presented in §[3] 



4.4. Simultaneous Effect of GC Mass, Color, and Size 

The best fits come from simultaneously fitting the effects of 
mass, color, and radius (Table |4] rows 6 and 8). Compared 
to just fitting the effects of mass and color, the fits are greatly 
improved {-lAip = 97.6), and the binned x^ values are rea- 
sonable (20.2 for both half-light radius and half-mass radius 
with 26 dof). Regardless of which size is chosen, the mass 
and radius exponents are the same. The addition of the radius 
actually increases the mass exponent beyond the fitting errors 
one derives if radius is ignored. The color index also differs 
when size is ignored. If the half-light radius is used, one finds 
a smaller index; however, we believe the correlation between 
half-light radius and color affects this fit, so that this does not 
give an accurate view of the dependence on the size of the 
cluster When the half-mass radius is used, one finds a larger 
index; however, the increase is only at a ^Icr level. The sizes 
of GCs play a critical role in not only affecting the number of 
LMXBs per GC, but also our interpretation of how strongly 
GC mass and color affect the number of LMXBs per GC. 

We believe that the simultaneous fitting of GC masses, col- 
ors, and half-mass radii provide the most physically motivated 
fit of the LMXB-GC connection. For this fit (Figure|9]l, we de- 
rive P = 1 .237!°;|j^^, S = 0.90^ \l, and e = -2.22!!^'. We note 
that by including the effect of GCs falsely matched to LMXBs 
we increased the mass exponent by 13%, increased the color 
index by 16% and decreased the radius exponent by 16%. 

If we require that the expected number of LMXBs scales 
linearly with the mass (/3 = 1), the negative log-likelihood 
function is increased by -2Aip = 10.7 for 1 more dof, and 
the binned x^ is 30.9 for 27 dof. Thus, we rule out a linear 
proportionality of mass at the 99.89% confidence limit. We 
discuss the implications of this in the next subsection. 

From Figure |9] we see that the declining probability for a 
GC with {g-z)> 1 .4 to contain an LMXB can be reproduced 
although we assume no break in the relation between color 
and A. This suggests that the masses and half-mass radii of 
GCs at these colors may be different than at bluer colors. In- 
deed, the fraction of GCs with r/, > 3 pc increases system- 
atically with color for (g-z) > 1.4. One possible explana- 
tion is contamination by diffuse ste llar clusters, which are 
redder, fainter co usins of GCs (e.g.. iLarsen & Brodiell2000t 
iPeng et aLll2006bh . Additionally, we note that the small varia- 
tions of the probability with galactocentric distance are repro- 
duced, even though the distance played no role in the fit. 

We used the indices from our best-fit to the Detected sample 
to determine the normalization for the Complete sample: 

/ . , \ 1.237 / ^ -2.22 

A, = 8.0 X 10-2 ( ] l0«"°(«-» ( ^ ) . 

(10) 

Since our method predicts the expected number of LMXBs 
per GC, we can determine the likelihood that a GC would 
contain multiple LMXBs. Note that M15 in our Galaxy has 
two active LMXBs (White & Angelini 2001); however, the 
LMXBs in M15 have Lx ^ 10^* ergs"' and would not be de- 
tected at the distances of the galaxies in our sample. We cal- 
culate the number of GCs containing multiple LMXBs, A'muit, 
as 

N^uU = J2^l-a + Xt)e-^'], (11) 

X.nX 

obtaining A^niuit = 1-3. Normalizing this to the number of 
GCs with LMXBs in the Complete sample (61), we esti- 
mate that 2.2+^-^% of GCs with LMXBs actually contain mul- 
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Fig. 8. — Identical to Figure |6] except with the results of fitting the dependence of the expected number of LMXBs on GC mass and color (Table|4] row 4) 
overlaid. Note that fitting the effects of mass and color alone does not reproduce the observed dependence of the probability on radius or color very well. 




Fig. 9. — Identical to Figure|6] except with the results of fitting the dependence of the expected number of LMXBs on GC mass, color, and half-mass radius 
(Table \4\ row 8) overlaid. Note that simultaneously fitting the effects of GC mass, color, and half-mass radius reproduces the observed dependence of the 
probability on all of the observed parameters reasonably well. 
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tiple LMXBs such that their combined luminosity is above 
3.2 X 10^'^ergs"'. This fraction should increase as the LMXB 
luminosity limit of the sample decreases. For example, the 
same calculation for our Detected sample with lower lumi- 
nosity LMXBs gives an expected fraction of S.2tfj% of GCs 
with multiple LMXBs. 

4.5. Implications of Simultaneous Fit 

For multivariate datasets, principal component analysis de- 
termines the set of orthogonal eigenvectors, linear combina- 
tions of the original variables, that minimize the Euclidea n 
distances to each eigenvector (e.g., iMurtagh & Heck|[l987h . 
The normalized variances of the data projected onto the eigen- 
vectors, the eigenvalues, are used to order the eigenvec- 
tors. The eigenvectors with larger eigenvalues are the best-fit 
eigenvectors to the data. Often, only eigenvectors with eigen- 
values greater than unity (i.e., those eigenvectors that are as 
good fits as the original variables) are retained. Principal com- 
ponent analysis among the masses, colors, and radii of GCs 
with LMXBs found two eigenvectors whose eigenvalues were 
above unity. The first eigenvector has an eigenvalue of 1.22 
and is dominated by an approximately equal combination of 
mass and half-mass radius, while the second eigenvector has 
an eigenvalue of 1.01 and is dominated by color. Therefore, 
we break our interpretation of the fits into two parts: an in- 
terpretation of the variation with color as a metallicity effect, 
and a discussion of the combined mass and half-mass radius 
variation as evidence for the dynamical formation of LMXBs 
in GCs. 

4.5.1. Metallicity Effect on the Presence of LMXBs 
We find a monotonic increase in the expected number of 
LMXBs with GC color of the form A> oc 10'' ''°-» "<«-<•). Us- 
ing data from NGC 4486, lACSVCSSl found that the prob- 
ability that a GC contains an LMXB is Px oc lO^-^^-S^^te-^). 
For small values of the expected number (A, ^1), these two 
quantities are nearly identical (Px ~ \, eqn. ||8l). Thus, our 
slightly different techniques are in go od agreemen t with re- 
spect to the fit to color variation. In lACSVCSS l the color 
variation was directly converted into a m etallicity variation 
by applying the lBruzual & Cha riot (2003) models, which give 
(£-z)~ (0.38 ±0.05) log and therefore Px ocZ^ ^^*" '. In 
this paper, we conver t ed colo r to metallicity following equa- 
tion (2) of iPeng et all (l2006ah . 

-6.21 + (5.14±0.67)(g-z) for^-z<1.05; 
^°^Zq " \ -2.75 + (1.83±0.23)(g-z) for§-z> 1.05. 

(12) 

We then replaced the 10^* term from equation |6] with 
{Z/Zq/ . We find 5' = 0.391!|]:;j^°, with no change to the 
variation with mass or radius. The d ependence on Z we find 
here is consistent with that found in IACSVCS3I and with re- 
cent results on the metallicity dependence of t he production 
of lik ely quiescent LMXBs in Galactic GCs (iHeinke et alJ 
120061) . Compared to fitting the effect of mass, radius, and 
color, the mass, radius, and metallicity fit is very slightly 
worse (-2A?/; = 2.0); however, the binned value is still rea- 
sonable (20.3 for half-mass radius, both with 26 dof). 

There are several ideas that could explain the relation be- 
tween metallicity and the likelihood that a G C will cont ain an 
LMXB, most of which were summarized in ACSVCS3L First, 
metal-rich stars may have larger radii and masses compared to 
metal-poor stars dBellazzini et al.|[l995b . which would make 



it easier to form LMXBs. This would increase both the num- 
ber of Roche-lobe overflow sys tems and the number o f bi- 
nary systems formed. Following lMaccarone et alj (l2004 l). we 
have estimated that the number of Roche-lobe overflow sys- 
tems has a weak metallicity dependence, Z" The number 
of NS binaries formed by tidal capture goes as Z" '^, but may 
have a lower exponent for exchange interactions. Thus, this 
model predicts a Z~''^ dependence and is inconsistent with 
o ur fit. 

Ilvanoval (l2006h suggested that the link between metallic- 
ity and outer convective zones was responsible for the in- 
creased formation efficiency of LMXBs in metal-rich GC. In 
this model, metal-poor main-sequence stars lack an outer con- 
vective zone. This absence turns off magnetic braking. Since 
magnetic braking leads to the strong orbital shrinkage that 
eventually leads to mass transfer as an LMXB, metal-poor 
stars fail to dynamically form many of the main-sequence/NS 
binaries that can appear as bright LMXBs. Although this 
model is intriguing, more work must be done to see if the 
removal of the outer convective zone would predict a depen- 
dence of the expected number of LMXBs on the GC metallic- 
ity that is consistent with our best-fit power-law dependence. 

If metal-rich GCs produce more NSs and BHs per unit 
mass, then the larger number of NSs and BHs could increase 
the number of LMXBs that can form. For instance, the initial 
mass function (I MF) could v ary with metallicity CGrindlavl 
11981) . Following IACSVCS3L we assume a power-law IMF, 
with the number of stars with masses between m and m + dm 
varying as N{m) cx m~''dm. We also assume that the IMF slope 
depends on metallicity, dx / d[Fe/H] =A. The IMF metallic- 
ity dependence must be A = -0.32^^■f■QJldex~^ to account for 
the m etallicity de pendence of the number of LMXBs per GC. 
As in ACSVCS3, the IMF dependence is in rough agreement 
with the analysis of LMXBs in the Galactic and M3 1 GC sys- 
tems, A -0.4 dex"' dBellazzini et al.lll995h . and IMF slope 
derived for Galactic GCs, A « -0.5dex~' (^Djorgovs ki et al.l 
Il993h . We conclude that a mildly metallicity dependent 
IMF can explain the metallicity dependence of the number of 
LMXBs per GC. However, we caution that there is evidence 
that massive star IMFs are metallicity indep endent (e.g., IMFs 
from OB associations: iMassev et al.iri995h . 

Irradiation-induced winds, which would be weaker in 
metal-rich stars due to more efficient metal line cooling, may 
affect the number of LMXBs observed in systems that are 
richer in metals. iMaccarone et alJ (l2004l) presents two ba- 
sic models that consider the effect of these winds. In the 
first model, the wind is ejected at the escape velocity of the 
donor star At the same accretion rate (i.e.. X-ray luminosity), 
stars that are richer in metals will have lower mass-loss and 
longer lifetimes. The rough predicted metallicity dependence 
for the number of LMXBs per GC is appro ximately Z"-^^, al- 
though we note that in making this estimate IMaccarone et al.l 
((2004) consider only the cooling function while the strength 
of stellar winds due to irra diation depends on both heating 
and cooling processes (e.g.. lBasko & Sunvaevlll9'73b . In the 
second model, a less dense, higher velocity wind is ejected. 
In this model, the mass-loss rates (i.e., lifetimes) of stars that 
are poorer in metals and stars that are richer in metals are 
about the same, but the accretion rates are not. Stars that are 
richer in metals are more luminous and the predicted metallic- 
ity dependence for the number of LMXBs per GC is roughly 
2~o.7(/3ir-i)^ where (3if is the differential X-ra y luminosity func- 
tion slope. Typical values of f3\f are 1.8-2.2 dKim & Fabbiand 
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Fig. 10. — Two-dimensional confidence intervals (at \a, 90%, and 2a) of the exponents of mass (f3) and radius (e) in the fit to the expected number of LMXBs 
per GC, XocMl^ 10*<s-=) r^^^^. The cross indicates the best-fit values {3 and e. The overlaid lines indicate the dependences of three likely dynamical-formation 

properties that depend on mass and radius: A; oc Fj^ , A, oc F;, „ = F;, p" , and A, oc . The power-law dependence on F;, is the most consistent, and a power-law 
dependence on S can be ruled out. 

form of the encounter rate'^). The mass and radius exponents 
clearly indicate that powers of the stellar crossing rate do not 

fit the data. 

The IACSVCS3I form of the encounter rate cannot be ruled 
out when considering the 2-dimensional confidence interval. 
However, fitting 



12004 for a single power law; however, the X-ray luminosity 
function may be consistent with a broken power law. At the 
higher luminosities (Lx > 5 x 10^'^ ergs"'), /?if = 2.8 ±0.6. At 
the lower luminosities that are characteristic of most LMXBs, 
f3\f =1.8 ±0.2. Thus a wide range of metallicity exponents 
are possible, but most are steeper than Z"''^ "'. Either of these 
models may be consistent with our fit; however, more detailed 
work developing irradiation-induced wind models is neces- 
sary to predict the metallicity dependence accurately. We note 
that systems with stronger irradiation-induced winds are pre- 
dicted to have more absorption and that spectral analysis by 
X-ray colors in Kim et al.. (2006) do not find evidence for this 
absorption. An ongoing similar analysis for the sample in this 
paper, as well as more detailed X-ray spectral fitting, should 
shed additional light on this issue. 

Finally, we note that many of the above models assume 
that the compact object in the LMXB binaries is a NS. 
However, our Complete sample contains GC-LMXBs above 
3.2 X 10-"^ ergs"'. Such sources are above the Eddington limit 
for a hydrogen-accreting neutron star, and the compact object 
in the system might be a BH. These models need to be devel- 
oped to include the effects of such binary systems. 



4.5.2. Dynamical Formation of LMXBs in GCs 

Since the efficiency of the dynamical formation of LMXBs 
in GCs is likely to depend on a combination of the GC mass 
and size, we plot the two-dimensional confidence interval of 
the exponents for both properties in Figure [TO] The effect of 
GC color was also included in these fits, but is not shown. On 
top of the figure, we have overlaid three lines: A, oc (powers 
of the stellar crossing rate). A, cx (powers of the encounter 



rate), and A, cx r,,,„ = r,,p" = M^^^"' r] 



-5/2-3a 
/7,cor 



(the lACSVCSH 



M 



3/2+a 



Ipc 



-5/2-3Q 



z 



(13) 



-0.200!!;;j^^ and 6' 



is a worse fit 



we find a 

(-2Aip = -4.4) than fitting mass, radius, and metallicity. With 
only one more dof, this fit is marginally (in)consistent at the 
96.5% confidence limit. The slightly worse fit is also sug- 
gested by the binned value of 33.56 for 27 dof. Any slight 
disagreement must be tempered by the presence of different 
assumptions in the analysis, namely the fact that IACSVCS3I 
did not assume a fixed c for all GCs as we do here. We do note 
that the value of a ~ -0.2 is roughly consistent with analysis 
of pulsars in Galactic GCs (Johnston et al. 1992). We display 
the 2-dimensional confidence interval of a and S' in the left 
panel of Figure [TT] 

The best fit comes from powers of the encounter rate. When 
we fit 

5' 



(14) 



we find ri = 0.818!JJ:!J^^ and S' = 0.389!{;:g^^ has only a sUghtly 
worse fit (-lAtp = -0.6) than fitting mass, radius, and metal- 
licity. Since it has one more dof, this is our preferred fit. Its 
binned value was 22.2 for 27 dof. Again, we used the in- 
dices from our best-fit to the Detected sample to determine 

a in this equation is related to that in equation (13) of I ACS VCS31 

"vcsa. by «vcs3 = 1.5 + a 
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Fig. 1 1 . — (Left:) Two-dimensional confidence intervals (at Icr, 90%, and 2cr) of the exponents of density and metallicity in the fit to the expected number of 
LMXBs per GC, A oc T;, (M/r^^ ^^|,)" (Z/Zq)* . The cross indicates the best-fit values of a and 5'. (Right:) The confidence intervals for the exponents of the 
encounter rate and metallicity in the fit to the expected number of LMXBs per GC, A oc Fj^ (Z/Z0)* . The cross indicates the best-fit values of 5' and r}. 



the normalization for the Complete sample: 

, 0.818 / ^ \ 0.389 

A, = 4.1 X 10"^ ' ' ' 



lOV 



z 



(15) 



The dependence of A, on encounter rate matches the 
Galactic values found for low-luminosity X-ra y sources by 
iJohnston & Verbunti d 19961) of a ~ 0.2 and by iPoolev etal] 
(l2003h of A (X r''-^'*^^"^^. Si nce low-luminosity X-ray sources 
repre se nt a m ixture (e.g., iBassa et aP l2004t iHeinke et al.l 
120031 I2OO6I; lPoolevetal.1 l2003l) of quiescent LMXBs 
(qLMXBs), cataclysmic variables (CVs), and chromospher- 
ically active binaries (CABs), such samples likely trace 
additional source formation and evolution processes when 
compared to a sa mple of active LMXBs . For example, 
iPoolev et ai] (120031) and lHeinke et al.l(l2006h show indications 
that qLMXBs are more consistent with a linear F dependence 
than fainter, hard X- ray sources likely t o include CVs and 
CABs. For qLMXBs. [Poolev et al.l (l2003h fin ds A oc ri "±°-5 



To compare to qLMXBs from iHeinke et aP (l2006h . we use 
the IACSVCS3I form of the encounter rate. Our a is defined 
su ch that it is 1 .5 less than the Alpha used in Figure 8, right 
of lHeinke et an r2006). and our 5' is their Delta. These values 
lie close to the > 50% probability contour, closer than a linear 
dependence on F without any metallicity dependence. How- 
ever, the > 50% probability contour value includes a large 
range of values both shallower (to the left of Alpha = 1.5) and 
steeper (to the right of Alpha = 1.5) than a linear dependence 
on F. The relatively limited numbers of Galactic qLMXBs in 
current samples make it difficult for these studies to constrain 
the dependence of the production rates of LMXBs in GCs. 

ACSVCS3 suggested two mechanisms to account for our 
observed less than linear dependence on the encounter rate; it 
might result from the competing rate of destruction of bina- 
ries, or it might be due to hardening (shrinking) of binaries. 

If both the formation and destruction rates of the close bina- 
ries per unit stellar mass) d epend on the encounter rate param- 
eter per unit mass F/,/M (Verbunt 2003) then in steady-state 
the encounter rate would cancel and the number of LMXBs 
per GC would just be proportional the number of stars, or 



A, cx M. This would be inconsistent with the dependencies 
of A, on mass, radius, and on F/, that we have found. How- 
ever, if the destruction rate of binaries due to encounters is 
comparable to but less than the destruction rate due to stellar 
evolution of the binaries, then an intermediate dependence on 
F/,/M (between linear and none) might result. Detailed simu- 
lations of binary formation, evolution, and destruction in GCs 
are needed to test this hypothesis. 

During exchange interactions, the orbits of wide binaries 
can be red uced, particularly in denser clusters. Following 
lACSVCS3l the binary encounter cross-section, C, can be writ- 
ten as a combination of the cross-section for tidal capture, 
a2, and the cross section for binary exchange interactions, 
o'3 C(p)^ where C,{p) encapsulates the effect of hardening as 
a function of density. If binary hardening is responsible for 
a = -0.2, then C = a2 + C(P)'^i oc p""-^. This would yield 
C(p) = (/i - cr2)/cr3, where /i is a constant, as a constraint 
on predictions for binary hardening. 

5. CONCLUSIONS 

We have compared the masses, colors, sizes, and galacto- 
centric distances of 270 GCs that contain LMXBs and 6488 
GCs that do not contain LMXBs in a sample of eleven galax- 
ies. There are clear differences between the masses, col- 
ors, and radii of GCs in these two samples, while the role 
of galactocentric distance has at most a weak effect. There 
is no significant difference between the masses, colors, or 
sizes of 61 GCs that contain LMXBs in the Complete sample 
(Lx > 3.2 X 10^^ ergs"') and 209 GCs with fainter LMXBs; 
thus, there is no evidence that luminous and fainter LMXBs 
are affected differently by the properties of the GCs they in- 
habit. 

We clearly show that GCs that are more massive and redder 
are more likely to contain LMXBs, confirming previous find- 
ings. Although we find that red-GCs are 3. 15 ±0.54 times 
more likely to have LMXBs than blue-GCs, we show that 
the detailed dependence changes continuously with GC color, 
rather than being a simple function of the overall population 
to which the GC belongs. 
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With the half-mass radius, we calcu late the relaxation 
timescale. In contrast to Galactic GCs (lB06l) . we find that 
a large number of GCs that contain LMXBs (~ 15%) have re- 
laxation times > 2.5 Gyr. Based on the results from our sam- 
ple, it does not appear necessary for GCs to survive for more 
than five relaxation timescales in order to produce LMXBs. 

Most notably, we find that GCs that have smaller half-light 
radii or half-mass radii are more likely to contain LMXBs. 
This paper presents the first clear indication that GC half-mass 
radius affects the likelihood a GC will contain an LMXB, 
although a size de pendence was implicit in the results of 
IJordan et al.ll2004bi 

Simultaneous fits of the dependence of the expected number 
of LMXBs per GC on the GC mass, color, and radius gave 



A, cxM' 



-0 075 



-z) -2-22!l 
/i.cor 



(16) 



Including the radius is important because fitting mass and 
color alone does not provide an adequate fit to the expected 
number of LMXBs, and the form of the dependence on mass 
and color changes when the effects of the size are included. 
The simplest model for LMXBs, that the number of LMXBs 
per GC is linearly proportionality to GC mass, can be ruled 
out at the 99.89% confidence limit. 
For our Complete sample, Lx > 3.2 x 10^'^ ergs"', we find 



A, = 8.0 X 10" 



M 



1.237 



10' 



,0.90 te- 



Ipc 



-2.22 



(17) 

We additionally provide for the first time an expression to esti- 
mate the expected number of multiple LMXB sources in GCs, 
which predicts that 2.2+j 7% of GCs with LMXBs actually 
contain multiple LMXBs such that their combined luminosity 
is above 3.2 x 10-"^ ergs"'. Thus, we predict that most GCs 
with high X-ray luminosities contain a single LMXB. 

Principle components analysis is used to show that the de- 
pendence of the expected number of LMXBs per GC on mass, 
color, and radius is mainly due to a dependence on a combina- 
tion of mass and radius, and a dependence on color We show 
that this result implies that the dependence on mass, color, and 
size is essentially equivalent to a dependence on the encounter 
rate F/, and the metallicity Z. The best-fit form is 



X,=A 



10^/ 



0.818: 



+0.050 
-0.049 



0.389: 



,+0.071 
0,067 



(18) 



This provides strong evidence for dynamical formation play- 
ing a primary role in forming LMXBs in the dense stel- 
lar environs of extragalactic GCs, strengthening the evidence 
presented in Jordan et al. (2004b). For our Complete sam- 
ple, Lx > 3.2 X 10^*^ ergs"', we find that the normaUzation is 
A = 4.1 X 10"2. 

The F/, exponent is consistent with that deri ved for low- 
luminosity X-ray sources in Galactic GCs ( Poolev et al.l 
|2003|), but inconsistent with the simple theoretical prediction 



that the number of LMXBs per GC be Unearly proportional to 
F/,. The hardening of binaries could explain the shallower 
dependence we observe. Alternatively, our use of F/, as a 
proxy for the encounter rate may affect the detailed depen- 
dence, particularly if core-collapsed extragalactic GCs prefer- 
entially contain LMXBs. Our ongoing study of the GCs and 
LMXBs in Centaurus A will test our use of F/,. 

The metallicity exponent is most consistent with either a 
metallicity dependent variation in the number of NSs/BHs per 
GC, such as on e might find in a metallicity dependent IMF 
jGrindlavll 19871). or effe cts from an irradiation induced wind 
(Maccarone et al.ir2004l) . although more detailed work devel- 
oping irradiation-induced wind models is necessary to predict 
the metallici ty depen dence accurately. The magnetic brak- 
ing model of llvanova! ("2006 ) may be consistent; however, this 
model needs to be developed over a range of metallicities. 

The combination of ongoing and upcoming deep CXO ob- 
servations and wide-field HST-ACS observations of NGC 
3379, 4278, 4365, 4697, and 5128 will be essential in extend- 
ing our ability to test the LMXB-GC connection to lower X- 
ray luminosities across entire galaxies. The greater number of 
GC-LMXBs will more tightly constrain the GC-LMXB con- 
nection and the theoretical models to explain it. The spatial 
distribution of GC-LMXBs and field-LMXBs may determine 
the level at which GC-LMXBs that have since been removed 
from their GC birthplaces contribute to the field population of 
LMXBs. It would also be very useful to extend the studies 
of LMXBs in early-type galaxies to include more lenticular 
galaxies. Our understanding of difference in the GC-LMXB 
connection between elliptical and lenticular galaxies is lim- 
ited by the small number of lenticular galaxies observed with 
CXO. This limits our ability to use LMXBs to trace the star- 
formation history of galaxies. A complementary approach is 
required that combines deep studies of individual lenticular 
galaxies, as in our approved study of NGC 1023, and shal- 
lower studies of larger samples of lenticular galaxies. 
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